In this paper, we suggest a simple iterative method to find the geodesic path on a surface parameterized by orthogonal curvilinear system between TWO given POINTs based on solving BOUNDARY VALUE PROBLEM. In this supposed method, an iterative algorithm is used for finding the sufficient initial VALUEs as the destination POINT agree with the BOUNDARY conditions. Geodesic determination between TWO given POINTs is formulated for a general surface, and specially tested for reference ellipsoid which has many applications in geosciences and geodesy. Accuracy of the method is independent on the distant between TWO POINTs on the surface. Moreover, it can be used in aviation and sailings for finding the shortest path between start and destination POINTs.